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Abstract 

Hyperion  is  a  hyperspectral  sensor  on  board 
NASA’s  EO  -1  satellite.  Its  spatial  resolution  is 
about  30  meters  with  a  swath  of  ~7  Km.  Though 
Hyperion  was  not  designed  for  ocean  studies,  its 
unique  spectral  configuration  (430  nm  -2400  nm 
with  a~10nm  step)  makes  it  especially  attractive  to 
study  the  effectiveness  of  such  kind  of  sensor  for 
observing  complex  coastal  waters.  In  this  study, 
Hyperion  data  over  two  sites  of  the  Florida  coasts 
were  acquired,  with  one  focused  on  the  clear  Key 
West  waters,  and  the  other  focused  on  the 
relatively  turbid  Tampa  Bay  waters.  From  both 
data  sets,  water  properties  and  bottom  bathymetry 
were  simultaneously  derived  from  atmosphere  - 
corrected  Hyperion  data  using  a  spectral  matching 
technique.  More  importantly,  in  the  top-to-bottom 
processing  of  Hyperion  data,  there  was  no  use  of 
any  a  prior  or  ground  truth  information.  For  the 
Key  West  site,  derived  bathymetry  and  water 
properties  were  validated  with  NAVOCEANO 
CHARTS  (active  bathymetric  LIDAR  system)  and 
field  measurements,  respectively.  It  is  found  that 
the  retrieved  depths  (in  a  range  of-  1  -  20  m) 
match  LIDAR  depths  very  well  (-15%  average 
error),  indicating  significant  potential  of  using 
hyperspectral  stellite  sensor  for  efficient  and 
repetitive  observation  of  shallow  coastal  regions. 

1.  Introduction 

For  many  Navy  operations,  environmental 
information  of  coastal  waters  is  critical  and  impacts 
the  success  or  failure  of  the  mission.  In  shallow 
waters  (less  than  20  meters)  the  optical  properties, 
bottom  reflectivity,  and  depth  are  all  essential 
elements  in  the  performance  of  mine  Electro- 
optical  Identification  (EOID)  systems.  The  Navy 
also  uses  active  electro  -optical  sensors  for  high 
resolution  surveys  of  coastal  bathymetry  and 
shoreline  characteristics.  Systems  such  as  the 
CHARTS  (Compact  Hydrographic  Airborne  Rapid 
Total  Survey)  sensors  employ  both  active  and 


passive  high  resolution  sensors  in  the  visible  range 
of  the  spectrum.  However,  knowing  when  and 
where  conditions  are  favorable  for  utilization  of 
such  assets  effects  cost  and  efficiency  of  all 
electro-optical  sensors. 

For  shallow  coastal  waters,  Lee  and  Carder  1  have 
demonstrated  that  reliable  derivation  of  water  and 
bottom  properties  from  spectral  remote  sensing 
requires  a  sensor  with  hyperspectral  capability. 
Current  operational  satellite  sensors  designed  for 
observation  of  water  properties,  such  as  SeaWiFS 
or  MODIS,  however,  have  few  spectral  bands 
along  with  a  large  spatial  footprint.  Such  sensors 
are  unable  to  provide  the  detailed  and  accurate 
information  for  shallow  coastal  environments. 

The  Flyperion  sensor  on  board  the  EO- 1  platform  2, 
designed  for  land  observations,  however,  has 
wavelength  bands  from  -430  nm  to  2400  nm  with 
a~10  nm  resolution.  Also,  the  ground  resolution  of 
the  sensor  is  about  30  meters.  Such  spectral  and 
spatial  characteristics  are  significantly  better  suited 
for  the  study  of  coastal  waters  than  the  sensors 
such  as  SeaWiFS  or  MODIS.  But  the  limitation  of 
this  sensor  is  its  low  signal-to-noise  ratios 
compared  to  SeaWiFS  or  MODIS. 

The  low  signal-to-noise  ratio  has  big  impact  on 
water  observation.  Water  targets  typically  have 
much  weaker  signals  than  land  targets.  Hyperion, 
being  designed  for  land  operations,  did  not  require 
a  high  signal-to-noise  ratio.  Hyperion  then  may  not 
have  enough  sensitivity  to  “differentiate”  the  subtle 
change  of  water  properties.  Consequently  it  was 
perceived  that  Hyperion  would  have  little 
usefulness  for  water  studies.  For  many  shallow 
coastal  area,  however,  due  to  the  increased 
turbidity  of  water  and  strong  reflectance  from  the 
bottom,  the  signals  emanating  from  the  water  could 
be  much  stronger  than  that  from  oceanic  waters. 
Therefore,  as  indicated  in  an  earlier  study  by 
Brando  and  Dekker 3,  Hyperion  imagery  could  be 
very  useful  for  coastal  areas  and  for  Navy 
operations. 
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In  this  study,  using  Hyperion  data  collected  over 
the  Florida  Keys  as  an  example,  we  present  an 
innovative  and  practical  method  for  the  derivation 
of  remote-sensing  reflectance  (Rrs)  for  such  an 
imperfect  sensor.  Hyperion  Rrs  is  further  applied 
to  a  newly  developed  hyperspectral  optimization 
processing  scheme  (HOPE)  for  the  derivation  of 
water  and  bottom  properties.  Unlike  traditional 
empirical  regressions  that  require  many 
assumptions  and  ground  truth  data  to  derive 
bathymetry  from  Rrs4'5,  HOPE  is  fully  automated 
and  free  of  the  requirement  for  ground  information. 
Hyperion  derived  water  and  bottom  properties  are 
further  compared  with  in-situ  measurements  and 
shown  to  be  in  excellent  agreement.  These  results 
indicate  that  with  the  innovative  method  of 
deriving  Rrs  and  the  advanced  optimization 
algorithm  for  water  properties,  Hyperion  imagery 
can  provide  vital  and  reliable  information  at  least 
for  shallow  coastal  environments. 

2.  Data 

LI  A  Hyperion  data  over  Looe  Key  (Florida) 
collected  on  October  26,  2002  was  provided  by  the 
USGS.  Since  our  focus  is  on  water  and  bottom 
properties,  only  information  from  430  -  925  nm 
was  used.  The  image  was  centered  at  24°42’39’ 

(N),  81°22’15”  (W).  Figure  1  presents  the 
collection  area  and  a  subset  of  the  image.  This 
subset  includes  clear  oceanic  waters,  shallow  sandy 
bottom  waters,  and  complex  out  flows  from  nearby 
land  with  mangroves. 

During  the  collection  of  Hyperion  data,  bathymetry 
of  the  study  area  was  surveyed  with  the  SHOALS 
system  (an  earlier  version  of  CHARTS)  -  an  active 
bathymetric  LIDAR  system.  In-situ  measurements 
of  remote-sensing  reflectance  and  water’s  optical 
properties  were  also  collected  at  six  sites  in  the 
area  (see  red  dots  of  Fig.  1),  with  Rrs  measured  by  a 
custom-made  spectroradiometer  and  water 
absorption  coefficients  measured  by  AG9 
(Wetlabs,  Inc.). 

3.  Methodology  to  derive  remote  sensing 
reflectance 

To  analytically  derive  water  and/or  bottom 
properties  from  any  satellite  data,  the  first  step  is  to 
get  high  quality  data  of  remote-sensing  reflectance 
(Rrs),  which  is  defined  as  the  ratio  of  water-leaving 
radiance  (L_w)  to  downwelling  irradiance  just 
above  the  surface  (Ed).  It  is  Rrs  that  solely  contains 
water  and/or  bottom  information. 


In  general,  the  radiance  measured  by  a  sensor  at 
any  altitude  (L_t)  can  be  expressed  as 

L_t  =  L_sky  +  t*L_w,  ( 1) 

with  L_sky  for  contributions  from  the  sky  and  L_w 
for  contributions  from  below  the  water  surface,  t  is 
the  transmittance  of  L_w  from  sea  surface  to 
sensor  altitude.  From  the  definition  of  Rrs,  one 
obtains 

Rrs  =  (L_t  -  L_sky)/(t*Ed).  (2) 

To  obtain  from  L_t,  L_sky,  t,  and  Ed  must  be 
known.  Traditionally,  standard  methods  calculate 
L_sky,  t,  and  Ed  based  on  models  of  radiative 
transfer  in  the  atmosphere  6,1 .  Since  L_w  in  general 
makes  up  about  10%  of  L_t,  such  methods  require 
high  accuracy  in  the  measured  L_t  (error  of  within 
1%  for  ocean  applications)  and  high  accuracy  of 
aerosol  models.  The  Hyperion  sensor,  however,  has 
an  accuracy  of  within  5%  in  the  measured  L_t 8. 

The  L_t  error  may  cause  a  50%  error  in  L_w  which 
in  turn  can  cascade  into  significant  uncertainties  in 
derived  water  properties.  For  such  a  poorly 
calibrated  satellite  or  aircraft  sensor,  an  innovative 
method  is  required. 

To  meet  such  a  goal,  we  developed  a  practical  and 
reliable  technique  to  derive  Rrs  from  Hyperion  data. 
This  technique  is  an  extension  of  the  cloud-shadow 
method  developed  by  Reinersman  et  al.  9.  In  that 
method 9,  L_sky  is  calculated  from  pair  of  adjacent 
pixels  that  is  under  cloud  shadow  and  under  the 
Sun,  with  t  and  Ed  calculated  from  radiance 
transfer  with  models  for  aerosols.  It  also  required 
the  sensor  to  be  well  calibrated  radiometrically  and 
spectrally. 

In  our  technique,  we  calculate  L_sky  in  a  similar 
but  simplified  fashion  of  that  of  Reinersman  et  al. 9, 
but  t  and  Ed  are  evaluated  quite  differently.  In  our 
technique,  t  and  Ed  are  not  explicitly  derived. 
Instead,  the  product  of  t  and  Ed  is  estimated  using 
the  reflected  radiance  from  the  top  of  the  cloud. 
Specifically,  for  the  pixel  under  the  Sun,  its 
radiance  (L_tl)  is  expressed  as: 

L_tl  =  L_sky  1  +  t*Ed*Rrs.  (3) 

For  the  pixel  under  the  cloud,  its  radiance  (Lt_2)  is 
L_t2  =  L_sky2  +  t*Ed_sky*Rrs.  (4) 

Fig.2a  shows  the  locations  of  a  pair  of  such  pixels 
with  the  corresponding  L_tl  and  L_t2  from 
Hyperion  in  Fig.2b. 

If  we  assume  that  L_skyl  =  L_sky2  =  L_sky,  then 
one  can  get 

L_sky  =  L_tl  -  (L_tl  -  L_t2)/(1-Ed_sky/Ed).  (5) 
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The  ratio  of  Ed_sky/Ed  is  estimated  using  Radtran 
10,  L  sky  is  then  easily  calculated  from  Eq.5,  in 
units  that  is  the  same  as  L_t  from  the  sensor,  either 
in  raw  counts  or  units  of  absolute  radiance.  Fig.2c 
presents  the  calculated  L_sky  of  this  study. 

To  calculate  Rrs,  values  of  t*Ed  are  still  needed. 

For  this  component,  we  use  the  radiance  from  the 
cloud  top  to  make  the  estimation.  As  above, 
radiance  from  a  cloud  top  can  beexpressed  as 

Lt_cloud  =  L_sky  +  t*Ed*Rrs_cloud.  (6) 
Fig.2d  shows  averaged  Lt_cloud  of  a  few  randomly 
collected  cloud-top  radiance.  Lastly,  Rrs  at  any 
pixel  can  be  calculated  as 

Rrs  =  (L_t-  L_sky)/(Lt_cloud  -  L_sky)*Rrs_cloud. 

(?) 

Fig.2e  presents  the  calculated  Rrs  for  the  pixel 
under  the  Sun  (Fig. 2b),  with  a  spectrally  constant 
Rrs_cloud  value  of  0. 1 59  (a  value  derived  based  on 
the  absolute  radiance  of  L_cloud  used  in  the 
calculation).  This  spectral  Rrs  shows  reasonable 
values  and  spectral  shapes  as  commonly  observed 
in  literature.  And,  the  Rrs  shows  smooth  spectral 
variation  in  the  430  -  800  nm  range,  as  compared 
to  the  slightly  noise  Rrs  obtained  by  Brando  and 
Arnold 3,  though  both  were  high  signal  coastal 
waters. 

One  obvious  advantage  of  this  new  approach  is  that 
L_t,  L_sky,  and  L_cloud  are  all  measured  by  the 
same  sensor  at  the  same  time;  therefore  the  derived 
Rrs  has  no  dependence  on  the  accuracy  of  L_t, 
because  sensor’s  response  function  (calibration 
factor)  is  canceled  out.  Also,  since  L_sky  is  the 
smaller  spectrum  of  the  image,  Rrs  calculated  by 
this  method  will  normally  be  positive  in  the  short 
wavelengths.  This  method,  however,  does  depend 
on  the  assumption  that  L  sky  is  nearly  uniform  for 
the  small  study  area. 

By  applying  Eq.7  to  the  entire  Hyperion  image,  Rrs 
is  calculated  for  each  pixel  within  the  scene.  Fig. 3 
compares  the  Rrs  from  Hyperion  to  Rrs  from  in  - 
situ  measurements  using  this  technique.  Excellent 
agreement  between  the  Hyperion  Rrs  and  the  in- 
situ  Rrs  is  demonstrated  in  four  out  of  the  six  in- 
situ  stations  (including  the  offshore  clear  water 
station),  in  both  spectral  shape  and  spectral  values. 
However,  for  two  stations  (St.2  and  St.5),  while  the 
spectral  shape  of  the  two  Rrs  are  quite  consistent, 
the  Hperion  Rrs  are  significantly  higher  than  in-situ 
Rrs.  The  reasons  for  this  mismatch  are  not  clear 
and  further  analysis  is  required. 

4.  Retrieval  of  environmental  properties 


To  derive  properties  of  the  water  column  and 
bottom  from  Rrs,  we  applied  the  optimization 
approach  developed  by  Lee  et  al.  "■ 12.  Briefly,  the 
approach  models  spectral  Rrs  as  a  function  of  five 
independent  variables  for  optically  shallow  waters, 

RM  =  F(aw(\),bM>P>G>X,B,H) 
RM  =  F(aM,bb„(*7),P,G,X,B,H) 

RM  =  F(awan),bbM,P,G,X,B,H) 

(8) 

Here  P  and  G  are  absorption  coefficients  of 
phytoplankton  and  gelbstoff  at  440  nm  respectively; 
X  is  the  backscattering  coefficient  of  suspended 
particles  at  440  nm,  B  the  bottom  reflectance  at 
550  nm,  and  H  the  bottom  depth.  To  derive  the  five 
unknowns,  a  spectral  optimization  scheme  with 
computer  processing  code  (HOPE)  has  been 
developed.  By  varying  the  values  of  the  five 
unknowns,  the  five  unknown  components  are 
considered  derived  when  the  modeled  Rrs 
spectrum  best  matches  the  Hyperion  spectrum. 

5.  Results  and  discussion 
The  derived  water  and  bottom  properties  from 
Hyperion  Rrs  were  compared  with  those  from  other 
independent  measurements.  Figure  4a  presents  the 
derived  image  of  water’s  total  absorption 
coefficient  at  440  nm  (a  parameter  for  water 
turbidity).  As  expected  for  a  coastal  system  Fig.4a 
shows  systematic  increase  of  a(440)  from  offshore 
to  inshore  in  a  pattern  parallel  to  the  coastal  line. 

In  addition,  Fig.4a  shows  the  outflow  of  in-land 
waters  into  the  coastal  system.  As  demonstrated, 
for  such  a  small  area  (~  7  km  by  20  km)  this  water 
optical  property  is  not  homogeneous  but  varies  by 
an  order  of  magnitude.  If  a  method  to  retrieve 
bottom  properties  relies  on  the  assumption  of 
homogeneous  water  properties,  then  it  will  have 
difficulties  for  such  an  area.  On  the  other  hand,  for 
pixels  that  are  in  parallel  to  the  coast,  we  see  a 
nearly  constant  optical  property  (though 
progressively  increase  to  coastal  line),  most  likely 
a  result  of  mixing  from  tidal  flow  and  alongshore 
currents. 

Fig.4b  compares  Hyperion  a(440)  (total  absorption 
coefficient  at  440  nm)  with  in-situ  a(440).  For  the 
five  stations  that  have  a(440)  measurements,  the 
stations  with  lower  a(440)  values  compared  very 
well  (average  difference  is  about  10%)  between 
Hyperion  results  and  in-situ  measurements.  For 
St.5  and  St.6,  where  a(440)  values  are  about  an 
order  of  magnitude  higher  than  that  of  other 
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stations,  the  average  difference  in  a(440)  is  about 
25%.  Because  both  St.5  and  St.6  are  very  shallow 
in  depth  this  larger  difference  may  be  due  to 
increased  bottom  reflectance  which  dominated  that 
from  scattering  of  the  water  column  (see  following 
section).  With  increased  contribution  from  the 
bottom  there  is  a  much  smaller  signal  contribution 
from  the  water  column  and  the  accuracy  of 
retrieving  water  properties  is  then  limited. 

However,  for  such  complex  region,  the  derived 
a(440)  value  is  still  quite  consistent  with  those 
from  in-situ  measurements. 

In  contrast  to  the  image  spatial  patterns  shown  by 
a(440),  the  depth  image  from  Hyperion  is  much 
different  (Fig.5a).  Instead  of  quite  uniform  pattern 
in  the  lower  portion  of  the  image  shown  by  a(440), 
the  bathymetry  i  mage  shows  distinctive  patchiness 
for  that  area,  with  a  shallow  sandy  bar  and  a  deeper 
channel  clearly  depicted.  Also  clearly  shown  is  the 
deeper  ship  channel  parallel  to  the  coast  (known  as 
Hawk  channel),  and  the  progressively  shallower 
depth  approaching  the  shore.  Compared  to  the 
conventional  bathymetry  chart  (see  Fig.  la),  the 
bathymetry  derived  from  Hyperion  shows  excellent 
consistency  with  the  chart,  but  Hyperion 
bathymetry  provides  better  resolution  and  much 
more  details.  Such  information  is  very  important 
for  Navy  applications  and  coastal  navigation. 

Fig.5b  compares  Hyperion  bathymetry  with  the 
bottom  depth  measured  at  the  five  stations.  For  the 
range  of  -2-13  meters,  both  depths  agree  with  each 
other  very  well  (less  than  10%  difference  on 
average).  Further,  Hyperion  bathymetry  is 
compared  with  the  bathymetry  derived  for  the 
SHOALS  system  (Fig. 5c).  For  over  70,000  points 
that  have  measurements  from  both  systems  and 
covering  a  range  of -1  -  20  meters,  the  bathymetry 
from  the  two  systems  show  good  agreement 
(average  difference  is  -15%).  Note  that  the  spatial 
resolution  of  SHOALS  is  about  1  meter,  while  the 
spatial  resolution  of  Hyperion  is  about  30  meter,  so 
binning  of  the  lidar  data  is  necessary  in  order  to 
compare  the  two.  Also,  Hyperion  has  an 
uncertainty  of -100  m  in  ground  registration.  With 
these  sources  of  errors  in  mind,  an  overall  error  of 
-15%  in  bathymetry  from  Hyperion  is  quite 
promising,  suggesting  significant  potentials  of 
using  Hyperion  for  observation  of  coastal  regions 
in  both  water  colurm  and  bottom  properties. 

6.  Conclusions 

We  have  demonstrated  that  with  our  innovative 
approach  of  deriving  Rrs  from  Hyperion,  excellent 


Rrs  spectrum  could  be  obtained  from  Hyperion 
data.  This  method  avoids  the  rigid  requirement  of 
accurate  radiometric  calibration  of  the  sensor,  and 
overcome  commonly  encountered  problem  of 
negative  Rrs  in  the  blue  when  using  standard 
atmosphere  correction  algorithm. 

Feed  Hyperion  Rrs  to  an  advanced  algorithm  to 
derive  properties  of  water  column  and  bottom,  we 
also  retrieved  surprisingly  impressive  results  for  a 
sensor  that  was  not  designed  for  water  studies.  This 
has  tremendous  implications  for  other  future 
sensors  and  existing  sensors.  The  implication  is 
that  such  techniques  can  extend  the  utility  of 
sensors  primarily  designed  for  land  applications  for 
some  coastal  applications.  While  these  results  are 
not  meant  to  suggest  land  satellites  can  be 
automatically  be  used  for  water  studies.  It  does 
suggest  that  information  can  be  obtained  from  such 
sensors  that  can  benefit  coastal  applications  for  the 
military  and  commercial  applications. 
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Figure  1.  (a)  Location  of  the  Hyperion  path,  (b)  A  subset  of  Hyperion  image  with  locations  of  the  six 
stations  with  in-situ  measurements. 


Figure  2.  Demonstration  of  “  Cloud -Shadow  Technique.  Panel  (a)  shows  selection  of  a  pair  of  pixels  under 
the  Sun  and  under  cloud;  (b)  L_t  of  the  two  pixels;  (c)  derived  L_sky  from  this  pair;  (d)  L_t  from  top  of  the 
clouds;  and  (e)  calculated  Rrs  for  the  pixel  under  the  Sun  (Fig. 2a). 


Figure  3.  Rrs  from  Hyperion  (green  line)  compared  with  Rrs  from  field  measurements  (red  line)  for  the  six 
stations. 
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Figure  5.  (a)  Bathymetry  image  derived  from  Hyperion;  (b)  bathymetry  compared  to  in-situ  value  collected 
from  shipboard  measurements. 
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Depth  from  Hyperion  [m] 


Figure  5c.  Comparison  of  Hyperion  bathymetry  and  SHOALS  Iidar-measured  bathymetry. 
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